Critical Lines and Massive Phases in Quantum Spin Ladders with Dimerization 
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We determine the existence of critical lines in dimerized quantum spin ladders in their phase 
diagram of coupling constants using the finite-size DMRG algorithm. We consider both staggered 
and columnar dimerization patterns, and antiferromagnetic and ferromagnetic inter-leg couplings. 
The existence of critical phases depends on the precise combination of these patterns. The nature 
of the massive phases separating the critical lines are characterized with generalized string order 
parameters that determine their valence bond solid (VBS) content. 
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I. INTRODUCTION 

The issue of the existence of massive or critical phases 
in quantum spin systems has motivated a great deal of 
study in strongly correlated systems since the seminal 
work on the Haldane conjecture This issue became a 
central problem when more complicated arrays of spins 
chains, known as quantum spin ladders, were discovered 
experimentally in cuprate materials that exhibit high-Tc 
superconductivity when they are appropriately doped , 
i- 

From a more fundamental viewpoint, the study of crit- 
ical and massive phases in quantum spin ladders offer the 
possibility of a testing ground for studying complicated 
quantum many-body effects, that in some instances un- 
derly the physics of unconventional phases of matter [J] . 
A variety of spin ladders with different number of legs 
have been synthesized based on cuprate materials, like 
SrCu20^, Sr2Cu^Oz etc. Q or in other family of 
compounds La4+4nCus+2nOi4+sn 07 ^"^d they present 
typically antiferromagnetic rung couplings among the 
legs of the spin ladders. It is also possible to find mate- 
rials with less conventional ferromagnetic rung couplings 
as well like in certain chemical compounds like PNNNO 
and PIMNO Q. In addition, the new tools to study 
strongly correlated systems based on optical lattices open 
the possibility of implementing a variety of quantum spin 
systems including ladders Q, [lOj . 

The phase diagram of quantum spin ladders with stag- 
gered dimerization was conjectured on the basis of ana- 
lytical non-perturbative methods [ll[ like the non-linear 
sigma model (NLSM) complemented with additional in- 
formation from the weak and strong coupling limits in 
the rung coupling constants and dimerization parame- 
ters. Later ^a series of different approximate analyti- 
cal studies [l^, have been favorable for the exis- 
tence of a critical line in the simplest case of a 2-leg spin 
ladder with staggered dimerization. Also, some prelimi- 
nary numerical methods with the Lanczos algorithm , 
have shown support for this fact for ladders with 
small size. Despite these several studies, a complete non- 
perturbative numerical analysis of these staggered low 
dimensional quantum spin systems have remained as an 



open problem. 

In this paper we study a 3-leg quantum spin ladder 
with different types of rung couplings (either antifer- 
romagnetic or ferromagnetic) and different dimerization 
patterns. These ladders are complex enough so as to 
serve as paradigmatic examples for testing the conjec- 
tured phase diagrams [ll|. In order to achieve a conclu- 
sive answer to the conjectured phase diagrams for this 
system, we resort to a non- per turbative numerical tool 
like the DMRG method [ill , [13 , , [H , in partic- 
ular we employ its finite-size version based on the sweep- 
ing procedure to improve the convergence of the iterative 
steps. The characterization of the different phases sep- 
arated by the critical lines is performed with a DMRG 
calculation of generalized string order parameters (SOP) 
[2l| that were introduced to distinguish between massive 
phases in dimerized spin chains [22l | even when the local- 
ized spins in the chain take on half-integer values. These 
SOPs are extensions of the originally non-local vacuum- 
expectation values introduced for integer spins and the 
like m, [13. 

The elucidation of the existence of a critical line in the 
phase diagram of quantum spin ladders with staggered 
dimerization is not straightforward when the ladders have 
end points as in the open boundary conditions geome- 
try demanded by the standard DMRG method. Thus, 
we have to resort to numerical analysis of the low-lying 
spectrum of excitations in order to extract the correct 
gap in the bulk of the system when studying the uni- 
versal properties of these ladders in the thermodynamic 
limit (length going to infinity) [2^ . 

The string order parameter is a theoretical construct 
in condensed matter that allows the characterization of 
massive phases with a VBS structure. It has never been 
measured experimentally. However, with the engineering 
of optical lattices it would be possible to address a direct 
measurement of this important quantity p^ . 

This paper is organized as follows: in Sect|ll] we in- 
troduce the two patterns of dimerization in 3-lcg lad- 
ders, one being columnar ([1]) with ferromagnetic rung 
couplings, while the other is alternating ^ with antifer- 
romagnetic rung couplings. In the former, we establish 
a helpful connection with the S* = | antiferromagnetic 
alternating chain in a certain strong coupling limit. We 
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recall several conjetures about the phase diagram of both 
3-leg Heisenberg models [ll| that motivate their study 
with DMRG in order to clarify them. In Sect lllll we per- 
form the DMRG calculations of the low-energy gaps in 
both models of 3-leg ladders as a function of the vari- 
ous coupling constants. With this information we can 
establish the phase diagram and thus we establish the 
validity of the conjectured diagram In addition, we 
can give a precise location of the critical lines and we 
find qualitative differences between these critical lines in 
each model. In Sect llVI we introduce generalized string 
order parameter to characterize the nature of the mas- 
sive phases separated by the critical lines found in the 
previous section. These string orders are measured with 
DMRG techniques and we show that they are a valuable 
tool for detecting VBS states in dimcrizcd quantum spin 
ladders with different patterns of dimcrization. SectfVl 
is devoted to conclusions. In appendix we study with 
DMRG the S = ^ alternating spin ladder and its gen- 
eralized string order parameters. This case appears as 
a limiting case in the study of the ferromagnetic 3-leg 
ladder with columnar dimerization and it is used as a 
guiding example to find the phase diagram. 



a) J(l+v) 




J(l-y) 
b) J(l+y) 




J(l-Y) 



FIG. 1: Pictorial representation of the Hamiltonians cor- 
responding to the completely antiferromagnetic (J' > 0) 
model ((2)1 with alternated staggering and b) the ferromagnetic 
(J' < 0) model (HI with columnar staggering. 



followed by a weaker one, that is, the bond alternation 
follows a columnar pattern. The explicit Hamiltonian of 
this arrangement is 



II. QUANTUM HAMILTONIANS FOR 3-LEG 
SPIN LADDERS WITH DIMERIZATION 

One of the main interests established in [ll| was the 
existence of an interplay between 3-leg ladders with a 
columnar dimerization pattern and ferromagnetic inter- 
leg, or rung, couplings on one side, and 3-leg ladders with 
alternating dimerization and antiferromagnetic rung cou- 
plings, on the other. The point was that both arrange- 
ments of 3-leg ladders should exhibit critical lines, while 
neither their precise location was known nor the nature of 
the massive phases they separated was determined. This 
is the open problem that we address here by means of the 
DMRG method. To this end, we start introducing both 
arrays of ladders since they are the candidates to exhibit 
critical lines. 

Thus, we are mainly interested in studying a possible 
connection between two different arrangements of spins 
i forming a 3-leg ladder which interact via Heisenberg 
terms. Both of them combine bonds of different strength 
parameterized by a constant 7 and two different types of 
coupling constants, J for the Heisenberg interaction be- 
tween spins along the legs, and J' for similar interactions 
between the rungs of the legs. Since the physics of the 
problem depends only on the ratio J'/ J, the constant J 
will from now on be given a fixed antiferromagnetic value 
J = 1. The other constant J' will have positive sign in 
one of the models and negative in the other one. In addi- 
tion to this difference, both models will also differ in the 
staggering pattern as we shall discuss below. 

We next describe the model corresponding to the re- 
gion with J' < 0, with the bond alternation pattern such 
that every one of the three legs begins with a strong bond 



L-l 

f=l,2,3 i=l 

L L ^ ' 

+ J'Y. • + J'Y. ^^(2) • S,(3), 

where L denotes the longitudinal length of the 3-lcg lad- 
der, Si(^) is a spin-i operator located at the site i of the 
£-th leg with < z < L, 7 is the dimcrization parameter 
that sets the relative strength of the bonds that interact 
with Heisenberg coupling constant J along the leg of the 
ladder and J' on the rungs of the legs. 

Since the coupling among legs is ferromagnetic, we 
know that for values \J'\ ^ 1 this model converges to an 
effective 5 = 3/2 staggered spin chain. Studies with the 
NLcrM method predicts [2^ for this chain the existence 
of three critical points in the interval 7 G [—1,1] placed 
at 7c = ±2/3 and 7c = 0. Numerical studies showed 
that in fact these points correspond to 7 = ±0.42 and 
7 = [131, [11]. Since the S" = 3/2 dimerized spin chain 
gives us valuable information about the expected massive 
phases of the ladder in the strong ferromagnetic regime, 
and since there exist some subtleties regarding the string 
order parameters that may lead to confusion, we have 
included a numerical study of this chain in appendix [XI 

For the completely antiferromagnetic regime J' > 0, 
we will use another staggering pattern that differs with 
respect to the previous one. In this case only the first and 
third leg begin with a strong bond, while the second one 
begins with a weak one, that is, the bond alternation is 
not columnar anymore but still follows a regular pattern 
that is called alternating. The Hamiltonian of this model 
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e=i,2,3 i=i 

L 



(2) 



+ J'J2S,{1)- S,(2) + J'J2S,{2) ■ S,{3), 



with the same conventions as before. A pictorial repre- 
sentation of this Hamiltonian together with the Hamil- 
tonian of the ferromagnetic model is shown in fig. [TJ 

With this arrangement, it is also possible to effectively 
and approximately map the model onto a NLcrM and 
again this formalism predicts a critical behaviour in the 
phase diagram of the couplings J'/ J vs. 7 How- 
ever, this behaviour is only reliable in the strong cou- 
pling limit J'/J ^ l.s More specifically, it predicts a 
critical curve running from the point (7 = 2/3, J' = 0) 
to (7 = 1, J' = 4/5) and another one which is the mirror 
reflection of the latter with respect to the J' axis. These 
predictions however shall be considered only as qualita- 
tive approximations of the real behaviour of the system. 
In this particular case, it is evident that the critical line 
must cut the 7 axis exactly at 7 = since that point cor- 
responds to two decoupled S' = 1/2 staggered Hcisenberg 
chains. However, this behaviour is missed by the NLtrM 
technique. On the other hand, this model has not appar- 
ent limits which can give us a hint on the phases that give 
raise. In section IIIII however, our DMRG computations 
will give strong evidence of their nature. 

Despite the differences in both models, i.e. different 
sign of J' and different staggering pattern, there are var- 
ious features that connect them. First of all and more 
important is that at least in the line 7 = both models 
are constinuosly related as we vary J' < to J' > 0. On 
the other hand, according to NLcrM, the first model is 
critical only in the region J' > while the second one 
has only critical lines in the complementary part J' < 0. 
This dual-like behaviour combined with the expectation 
that the ground state of both models is a valence bond 
solid, rouse the belief about the possibility of establish- 
ing a connection among their phase diagrams [ll| . We 
shall see to what extend these expectations are fulfilled 
with the help of the DMRG technique and the general- 
ized string order parameters. 



III. GROUND STATE DEGENERACY AND 
EXISTENCE OF CRITICAL LINES 

Massive quantum phases are characterized by an en- 
ergy gap from the ground state (degenerate or not) to 
the first excited state. On the contrary critical phases 
are characterized by a gapless spectrum between these 
energy levels. We have used the finite-size DMRG al- 
gorithm to compute the low energy levels and thus the 
corresponding gaps in order to identify the gap in the 
bulk of the system when we send the length of the 3-leg 
ladders to infinity (thermodynamic limit). 



It is known that in the case of integer spin chains, some 
configurations of VBS states can break a hidden Z2 x Z2 
symmetry [2l| that makes the ground state degenerate. 
This degeneration is a reflection of the spin-end effects in 
a antiferromagnetic Heisenberg chain of S* = 1 spins [2^ . 
In systems other than the parent Hamiltonians of the 
VBS states, but close enough to this picture and when 
using open boundary conditions, typically this degener- 
acy is approximate and the energy of the near-degenerate 
states decays with the size exponentially to a unique in- 
finite volume ground state. 




FIG. 2: (Color online) Energy gaps between the ground state 
and two first excited states in model ((2)1 ( Up) and model 
(Down). For the sake of clarity A20 and A21 are shown only 
for one value of J'. Namely J' = 1.0 (Up) and J' = -0.4 
(Down). It can be seen how the gap A21 indeed represents 
properly the gap of the spectrum irrespective of the degener- 
ation of the ground state. 

Fig. [2]shows the energy differences between the ground 
state and two first excited states in both models ^ and 
We observe in this figure that Aio accounts for the 
degeneracy mentioned in the previous paragraph: in the 
completely antiferromagnetic model ([2|) the first excited 
state is clearly above the ground state in the phase corre- 
sponding to low values of 7, while it is degenerate in the 
rest of the 7 interval. On the contrary, the ground state 
of the ferromagnetic model ((T|) is degenerate for low val- 
ues of 7 while it has a finite gap in the rest. As for A20, 
it corresponds in fact to the gap of the spectrum in the 
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FIG. 3: Critical curves of model ([1} with J' > and al- 
ternated staggering, and model ([2} with J' < and colum- 
nar staggering. Each point in the critical lines has been ob- 
tained keeping one parameter fixed and finding the value of 
the other parameter that minimizes the gap A21. Compu- 
tations have been performed on ladders of size L = 3 x 150 
retaining m — 400( J' > region) and m = 450( J' < region) 
states of the density matrix. 



degenerate regions while it clearly differs from the gap 
Aio in the non degenerate regions. Finally, the energy 
difference A21 coincides in both regions and both models 
with the gap of the massive phases. This holds true up 
to slight deviations due to finite size effects and irrespec- 
tive of the degeneracy of the ground state, Thus, we may 
conclude that this gap A21 is in fact the same one that 
survives when using periodic boundary conditions. 

As mentioned in the previous section, the existence 
of critical lines (characterized by a gapless spectrum) in 
both models ([T]) and ([2]) is supported by arguments com- 
ing from the strong coupling limit in the case of model 
([T|), and NLcrM valid in both of them. These arguments 
are however not conclusive. In this section we will prove 
numerically that these lines exist and we give an accu- 
rate estimation of its shape and location. To this end, 
we will find the critical values "fc{J') that make the gap 
A2i(7c(^')) vanish. It is very important at this point 
to emphasize that an exactly vanishing value of the gap 
shall only be attained in the thermodynamic limit. For 
finite size systems, the magnitude of the gap remains fi- 
nite and gets closer and closer to zero as we increase the 
size. In our case however, critical points are separating 
massive gapped phases and therefore, for large enough 
sizes but still computationally feasible, the gap at these 
points attain a local minimum value and can be accu- 
rately computed. 

Fig[3]shows the critical region computed for both mod- 
els ([2|) and ([T]). In order to compute this curves, we have 
used the finite DMRG algorithm in ladders of L = 3 x 150 
sites. For the completely antiferromagnetic model ^ we 
retained m = 400 states of the density matrix and a Lanc- 
zos tolerance equal to I0~^. The ferromagnetic model ([T|) 
turned out to be numerically more dcmmanding and we 
set m = 450 and the tolerance equal to 10^^". Two 



sweeps of DMRG were enough in both models to make 
the energies converge. These results clearly confirm the 
conjectured phase diagrams for these 3-leg models [Tlj . 

The solid lines in fig. [3] are only a guide for the eye. 
We have however used our numerical data to stimate the 
best fit to that critical lines. For the region J' > we 
have used a simple potential function of the form 



j: 



(3) 



The best value of each parameter has been obtained 
performing a least square fit and are equal to a = 1.59 ± 
0.01 and r = 0.72 ± 0.01. As for the critical hne in the 
semiplane J' < 0, we have used for the region close to 
the vertical asymptota a relation of the form 



J' 



C 



iic - ay 



(4) 



And the values that best fit the data have been found 
to be C = 0.38 ± 0.02, s = 1.07 ± 0.02 and a = 0.427 ± 
0.001. Notice that the value of this last parameter is 
in good agreement with previous computations of the 
critical point of the S — 3/2 alternating dimcrized chain. 



IV. 



MASSIVE PHASES AND GENERALIZED 
STRING ORDER PARAMETERS 



(2,1)- VBS 



i ) (1,2)-VBS 



(3,0)- VBS 



FIG. 4: Valence bond solid diagrams of the phases that give 
raise the models discussed in this paper. Each small solid 
circle and line represents both spin- 1/2 variable and a singlet 
pair repectively. The large open circles represent the sym- 
metrization of the spin 1/2- variables on each leg to create a 
spin-3/2 variable. 

In this section we will characterize the quantum phases 
that appear in the phase diagram in FiglSl To achieve 
this goal, we will resort to the generalized string order 
parameter 21 1, [2^, which are able to detect the VBS 
state character of dimerized spin systems even when the 
local spins take on half-integer values. These parameters 
arc generalizations for arbitrary complex phase of the 
original string order [2^, [2^ parameter first proposed 
for the case of integer spin S = 1. Resorting to the 
VBS picture, massive phases corresponding to valence 
bond solids can be denoted according to the number of 
valence bonds formed with the contiguous sites, i.e, one 
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particular valence bond solid can be denoted as (to, n)- 
VBS with m + n = 2S. 

For instance, we have already mentioned in the previ- 
ous section that in the strong coupling limit, the colum- 
nar dimerized 3-leg ladder ^ effectively becomes a. S ~ 
alternating spin chain. Thus, in this case we have 
3. 

The definition of the generalized string order parame- 
ter extended to our particular three leg ladder with arbi- 
trary size L = 3 X £ is 



3 
2 

m + n 



2j-l 

hm {S^,cxp{ie ^ S^,)SI/, 

k=2i 



(5) 



with 0<i<j <e/2 and S'f = S'f (1) + Sf (2) -I- S'f (3). It 
is actually not necessary to consider j — i 3> 1 to obtain 
accurate values of the parameter and typically a value of 
i~ j of some few tens is enough to give values very close 
to the infinite limit value, considered that i and j are well 
within the bulk and far away from the edges. Thereby, 
for convenience we will work with the parameter defined 
as 

2i-i 

0,tr(2z, 2j, e, 9) = {S^.cM^O ^ Sl)SI^) (6) 

k=2i 

with i.j and defined as before. 

It has been shown [2l| that the generalized string order 
parameter evaluated in 9 = n acts as an order parame- 
ter since it vanishes or not depending on the number of 
bonds n being odd or even. Moreover, the shape of the 
string order parameter in the region 9 e [0, 27r] provides 
us with valuable information about the VBS character of 
the phases since the number of zeros in this range coin- 
cides with the number of bonds to [21| . 

Fig. [5] shows the parameter Ostr computed in the com- 
pletely antiferromagnetic model ^ in the whole range 
of 7 for various values of J' . The operator clearly dis- 
tinguishes regions where it is finite from others where it 
vanishes. Moreover for a fixed value of J', the value of 
7 where it decays to zero coincides [s^l with the critical 
value 7c corresponding to that value of J' in the critical 
line of Fig. [S] On the other hand. Fig. [6] shows Ostr 
computed in the ferromagnetic model ^ in the strong 
coupling regime. For strong values of the ferromagnetic 
coupling J' we should expect that our ladder behave like 
an effective S = 3/2 alternating spin chain. Indeed, we 
can observe that the string order parameter is clearly 
non-vanishing above 7 = 0.42, the critical point of the 
chain. As for the region below this point, the tendency 
of the string order parameter is to decay to zero as we 
increase the size of the system, except for the point 7 = 
and its vicinity. In fact, this behaviour is anomalous since 
7 = is critical in the chain and therefore the string or- 
der parameter should vanish. In the appendix we have 
addressed this issue with the pure 5 = 3/2 alternating 
dimerized chain. Our study in the chain explains the be- 
haviour of the ladder and shows that indeed the string 




FIG. 5: Generalized string order parameter Ostr(i = 20, j = 
42, £ — 80, 8 — Tv) computed in the completely antiferromag- 
netic model ([2]) with alternated staggering for several fixed 
values of J'. Regions where the string order parameter van- 
ishes correspond to a different quantum phase from that where 
it is non null. Notice that there exists a certain J' above which 
the system only exhibit one quantum phase irrespective of the 
value of 7. This value corresponds to the J' coordinate with 
7c = 1 of the critical line of figure |3l 



0.25 




FIG. 6: Generalized string order parameter Ostr(i = 1/2 — 
24, j = 1/2 -I- 26, £ = n,9 = -k) computed in the strong fer- 
romagnetic regime J' = —25 of the ferromagnetic model ([TJ 
with columnar staggering. See the text for explanations. 



order parameter decays to zero also at this point. The 
decay rate is however slower and it is not enough to in- 
crease the size of the ladder. Fig. [8] (down) shows that 
we have additionaly to consider sites i and j further and 
further apart to make the string order parameter decay 
at 7 = 0. 

With the considerations explained in the previous 
paragraph, we have proved that the string order param- 
eter ([6]) evaluated in = tt is valid to identify quantum 
phase transitions in both models ^ and Now we 
want to extend the study to the whole 9 domain to test 
the nature of the massive phases. In fact, from definition 
® Ostr{i,jAiO) = Osti{i,j,i,9 + Tr) and hence we can 
restrict the study to the range 9 G [0,7r]. We will study 
first the ferromagnetic model ([T|) with columnar stag- 
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FIG. 7: (Color online) Generalized string order parameter 
Ostr(i = 20, j = 42, £ = 80,0) computed in the completely 
antiferromagnetic model at 7 = 0.5 ( Up) and the ferromag- 
netic model at 7 = 0.35 (Down). According to fig. [3] this 
values cut the critical lines at (7c = 0.5, J'c — 0.96) and 
{7c = 0.35, J'c = —6.1) respectively. It can be observed that 
both graphs show a change in the number of zeros near these 
points. 



gering. As we commented in the previous sections, we 
can guess the phases of the diagram going to the strong 
coupling limit |J'| ^ 1. In this regime the ferromag- 
netic coupUng among rungs is the leading interaction and 
the ladder transforms into an effective 5 = 3/2 alternat- 
ing spin chain. Resorting to continuity arguments, the 
phases of the ladder must be the same that appears in 
the strong coupling limit, i.e, a (2, 1)-VBS in the region 
< 7 < 7c and a (3,0)-VBS when 7c < 7 < 1. In figEI 
(down) we show Ostr computed in the 9 domain. All 
the curves appearing in the figure correspond to a fixed 
7 = 0.35. According to Fig. [3] this value of 7 cuts the 
critical line in a certain value of J' and so we should no- 
tice a qualitative change in the curves in Fig. [7l This 
change can be observed since plots corresponding to very 
negative J' have a local minimum ai = n, while it 
changes to become a maximum as we move towards J' 
close to zero. Hence, the number of zeroes in the domain 
6 e [0, 27r) moves from one to two. In fact, the string 
order parameter is not strictly equal zero in our graphs, 
but this fact has been already pointed out in spin chains 



[28l | and conclusively proved that it was due to finite size 
effects. The values closer to zero are attained consider- 
ing larger sizes. According then to the VBS notation and 
our DMRG results, we can label the phases of model ^ 
as (2,1)-VBS in the region < 7 < 7^ and (3,0)-VBS 
for 7c < 7 < 1 as expected from the knowledge of the 
5 = 3/2 chain. 

The completely antiferromagnetic model is more diffi- 
cult to guess a priori the quantum phases that gives raise 
to or even if its ground states are valence bond solids. 
We have used again the generalized string order parame- 
ters to check the nature of the phases. FigjT] (up) shows 
the string order parameter as a function of 9 and a fixed 
value 7 = 0.5 which cuts the critical line. It can be seen 
that indeed the SOP behaves as expected for a valence 
bond solid VBS state and two phases can be identified 
attending to the number of zeros in this domain. For 
higher values of J' the SOP has a maximum at 6' = tt 
and only one zero in the region < 9 < 2tt. As we 
consider lower values of J', more precisely in the inter- 
val from J' = 1.0 to J' = 0.8, the SOP at 6* = tt falls 
abruptly to zero and therefore the SOP has two vanishing 
values in the aforementioned interval. From these results 
we can conclude firstly, that both massive phases can be 
properly described as valence bond solids and also they 
can be identified as a (1, 2)-VBS for < 7 < 7c and a 
(2, 1)-VBS for 7c < 7 < 1. 



V. CONCLUSIONS 

We have given a precise meaning to the conjectured 
phase diagrams [ll[ corresponding to 3-lcg Heisenberg 
ladders with columnar dimerization and ferromagnetic 
rung couplins on one side, and similarly for alternating 
dimerization with antiferromagnetic couplings among the 
rungs. Although both models exhibit critical lines, their 
qualitative form is different: in the former, the critical 
line approaches an asymptota at the critical value 7c = 
0.42 were it effectively becomes an alternating 5=1 
spin ladder. On the contrary, the latter model does not 
exhibit any asymptota but the critical line meets the wall 
7 = 1 of the phase diagram. 

Moreover, we have also clarified the valence-bond-solid 
nature of the massive phases separated by the critical 
lines in the phase diagram. In this regard, we have found 
that the generalized string order parameters are a very 
good tool for characterizing VBS state phases in lad- 
ders with a variety of dimerization patterns. Our results 
are based on extensive calculations using the finite-size 
DMRG technique. 
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131, L = 180,6* = tt) (solid line) and 0,7,™ = 52, j = 
130, L = 180,^ — 7r)(daslied line) computed in a S = 3/2 
alternating dimerized chain. Down: string order parame- 
ter Osti'^ at 7 = varying the total length of the chain 
OstT"(i = n/2-24, j = n/2+26, L = n,e^7r) (empty circles) 
and varying the distance i — j, Os^™(i = (L — n+l)/2 + l, j = 
i + n-l,L = 200, e = Tr) (solid circles). 
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APPENDIX A: THE S' = 3/2 STAGGERED SPIN 
CHAIN 

It was already pointed out by Oshikawa [2l[ that half- 
integer systems are amenable to have many order pa- 
rameters. We want to compare in this section two order 
parameters which despite differing slightly in their defini- 
tion have in fact quite different behaviour and may even 
lead to confusion. We will see that both parameters cap- 
ture the quantum phase transition of the chain, but only 
one can go further and give evidences of the quantum 
phases themselves. 

The string order parameters are defined as 



ot 



,L)= .hm (5|,exp(z0 5^ 51)^2^) (Al) 



l=1i 




FIG. 9: (Color online) string order parameter 0|^r''"(i = 
22, j = 80, L = 100, &) computed below and above the critical 
point 7c = 0.42 of a S = 3/2 alternating dimerized chain. 
The number of zeros of this parameter determines the nature 
of the valence bond solid at each phase. 



dd/ 



,£) 



lim 



2j 

E 

=2i+l 



2j+li 



with <i < j < L/2, and 



(A2) 

with <i < j < L/2. 

Notice that, including the edges, both definitions in- 
volve an odd number of spins between sites 2i and 2j or 
2i + 1 and 2j + 1 and an even number of antisymmetric 
operators under spin flip. This condition is required to 
obtain a value of the mean value different from zero. 

It has been commented in the previous sections that 
the S = 3/2 alternating dimerized chain has three critical 
points at 7c = and 7c = ±0.42 in the interval 7 £ 
[— 1 , 1] . The phase diagram is symmetric respect 7 
and therefore we can constrain our study to the region 
< 7. 

As regards the interval < 7 < 7c the ground state 
is known to be a (2, 1)-VBS while for 7c < 7 < 1 it is a 
(3,0)-VBS. The first issue we have to check is wether or 
not the parameters defined above can make explicit this 
quantum phase transition. In fig. [5] (up) we have plotted 
both parameters 0°^/^{tt, L) and 0°^°"(7r, L) in the range 
< 7 < 1. This graph shows that both operators are 
finite at one side of the critical point while thay vanish at 
the other. We can conclude then that they act as proper 
order parameters in the phase transition. However a clear 
major difference can be noticed from this figure since 
each one of these order parameters in fact vasnishes on 
differents sides of the critical point. 

On the other hand, the behaviour of the string order 
parameters at the critical point 7 = 6 requires further 
insight. We know that critical points correspond to crit- 
ical ground states where the hidden order measured by 
string order parameters must vanish. In fig. [5] (up) it 
is not however clear that the value of the string order 
parameters decays to zero at this critical point. We have 
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FIG. 10: (Color online) string order parameter 0°ti'^{i = 
21, j = 79, L = 100,61) computed below (Up) and above 
(Down) the critical point 7c — 0.42 of a S = 3/2 alternat- 
ing dimerized chain. In this case the string order parameter 
in the 9 does not clearly determine the valence bond phases 
of the chain. Notice also the appreciable change of scale from 
one phase to the other. 



addressed deeper this issue in fig. [8] (down). This graph 
shows a finite size scaling analysis of the string order 
parameter 0°™" at the critical point 7 = 0. It can be 
observed that the scaUng and decay of the parameter is 
more influenced by the distance i— j than the total length 
of the chain. The rise of the string order parameter for 
large distances i—j close to the total lenght L is ex- 
pected due to finite size effects. For values far enough 
of the edges the tendency of the string order parameter 
is however to vanish as expected increasing the distance 



A regards the shape of these parameters for arbitrary 
6, fig. [S] shows 0°^^'^{9,L) for two different values of 7 
below and above the critical point. From the knowledge 
that we have of the phases of the chain we can see that 
this order parameter behaves as pointed out in pl| and 
the number of zeros identifies the valence bond phase 
ocurring. In effect, for values 7 < 7c the number of 
vanishing values in < < 27r is two and coincides with 
a (2, 1)-VBS while for 7 > 7c there are three null values 
corresponding to the (3, 0)-VBS phase. 



On the other hand 0°^°'^{9,L) is plotted for various 
values of 7 in fig. [TOl Two features can be remarked 
of this parameter: it also shows a qualitative change in 
its shape as wc move from one phase to the other. In 
this case however the relation of the number of zeros 
with the nature of the phase is not clear. Besides, the 
characteristic scale of the parameter differs significantly 
in both phases. 
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